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Abstract 



We study neutrino masses in the framework of the supersymmetric inverse see- 
saw model. Different from the non-supersymmetric version a minimal realization 
with just one pair of singlets is sufficient to explain all neutrino data. We com- 
pute the neutrino mass matrix up to 1-loop order and show how neutrino data 
can be described in terms of the model parameters. We then calculate rates for 
lepton flavour violating (LFV) processes, such as u — > &y, and chargino decays 
to singlet scalar neutrinos. The latter decays are potentially observable at the 
LHC and show a characteristic decay pattern dictated by the same parameters 
which generate the observed large neutrino angles. 
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1 Introduction 



Currently there are only very few indications for physics beyond the standard model 
(SM), the most important ones coming from neutrino physics and cosmology On the 
one hand, neutrino oscillation experiments [1] have shown that at least two neutrinos 
have non-zero masses and that mixing angles in the lepton sector are surprisingly large 
[2]. On the other hand, data from the WMAP satellite [31 H] and large scale structure 
formation [5] have provided convincing evidence for the existence of non-baryonic dark 
matter. 

In this paper we study a minimal supersymmetric version of the inverse seesaw 
[6]. The model is capable to explain all neutrino data with only one pair of singlet 
superfields. It contains a new dark matter candidate - the scalar singlet - which can 
give the correct relic density [7] and it gives potentially testable predictions for both, 
supersymmetric phenomenology at the LHC and low energy lepton flavour violating 
decays, such as /i — > ej. 

Neutrino masses are not part of the SM, but models which can explain oscillation 
data are quite easily constructed. Indeed, it was pointed out already in [8] that for 
Majorana neutrinos the mass matrix is described by a unique dimension-5 operator, 

m u = L { HL){HL). (1) 

All models which reduce to the SM particle content at low energy are merely different 
realizations of this operator and at tree-level there are just three basic contractions 
which give rise to eq. ([1]) [9]. The literature is completely dominated by only one of 
them, based on the exchange of heavy singlets [TUj [TTj . This is now commonly called 
the (type-I) seesaw mechanism. 

In type-I seesaw the smallness of the observed neutrino masses is attributed to 
the large mass of the singlets (z/ c ) and for / ~ 0(1) current neutrino data indicates 
A ~ 10 15 GeV. Obviously, if this ansatz is the correct explanation for neutrino masses, 
it will never be directly testecj^. However, the smallness of m u could be understood 



^^Dirac neutrinos can just as easily explain oscillation data. However, Dirac neutrinos require 
Yukawa couplings of order 0(1O -12 ) or smaller, thus there is no conceivable experimental phenomenol- 
ogy outside the neutrino sector for Dirac neutrinos either. 
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as well, if / is / <C 1. The classical examples for this situation are radiative neutrino 
mass models [121 [TBI [Hj . 

In the inverse seesaw model [6] the particle content of the SM is extended by one 
or more pairs of singlets, call them v c and S, which form "heavy" pseudo-Dirac pairs. 
The smallness of m v is then attributed to a small lepton number breaking parameter, 
fig. The smallness of this parameter is natural in the t'Hooft sense [15] . since a 
vanishing fig restores a symmetry of the theory. Similar to the ordinary type-I seesaw, 
in the inverse seesaw only one non-zero neutrino mass for the light neutrino fields is 
generated for each pair of singlets. A non-supersymmetric inverse seesaw thus needs at 
least two pairs of singlets to explain neutrino oscillation data [16] . As we show below, 
in a supersymmetric inverse seesaw one pair of singlets is sufficient to explain the 
experimental data. In such a minimal supersymmetric inverse seesaw model (MSISM) 
one neutrino mass is generated at tree-level, while a second non-zero mass is due to 
the scalar neutrino-antisneutrino loop [17]. The scheme we consider is reminiscent of 
bilinear R-parity violation, which is also of the hybrid "tree + loop" type [18] . 

Supersymmetrizing the inverse seesaw offers additional advantages^. Cosmology 
requires the existence of a non-baryonic dark matter (DM) candidate and SUSY with 
conserved R-parity offers a WIMP candidate in the form of the lightest supersymmetric 
particle (LSP), for reviews see for example [19} [20] . In the minimal supersymmetric 
extension of the standard model (MSSM) only the lightest neutralino remains as a 
CDM candidate, since left sneutrinos have been ruled out as cold dark matter by a 
combination of experimental data from LEP and direct detection experiments [21j . 
Right sneutrinos could be the CDM, however, for in the case of pure Dirac neutrinos 
as well as in the case of the standard type-I seesaw Majorana neutrinos, the sneutrinos 
are expected to have such small couplings to all ordinary particles that they can not be 
thermally produced dark matter. Non-thermal right sneutrino DM has been discussed 
in [221 [23] . Right sneutrinos could be thermalized in the early universe, if they have (a) 
enlarged left-right mixing [Ml EE); (b) a large quartic coupling to the Higgs fields [26] ; 
(c) an extra U'(l) under which sneutrinos are charged [27J or (d) within the NMSSM, if 
the sneutrinos have a large coupling to the NMSSM singlet [28]. In the supersymmetric 

2 A supersymmetric extension of the SM which adds only singlets inherits all the standard argu- 
ments in favour of SUSY, such as providing a (technical) solution to the gauge hierarchy problem, 
gauge coupling unification, etc. 
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inverse seesaw, which we consider here, the singlet scalars are expected to be thermal 
cold dark matter candidates [TJ, since the neutrino Yukawa couplings are much larger 
than in the standard type-I seesaw. 

The large Dirac neutrino couplings lead necessarily also to non-zero lepton flavour 
violating processes, such as fi — > cy and LFV supersymmetric particle decays. We 
therefore calculate BR(£j — > li + 7) and compare the expected rates with experimental 
sensitivities. If SUSY particles are light enough to be produced at the LHC, the new 
singlet states can appear in the decay chains, potentially altering the phenomenology. 
This is especially important in case one of the singlets is the LSP. We therefore also 
calculate the decays xt ~ * h + N a , where N a stands for a scalar neutrino. The flavour 
of the lepton in these decays can be tagged and traces the lepton flavour violating 
couplings of the sneutrinos. We show how these LFV couplings are related to the 
observed neutrino angles in the theoretically preferred part of the parameter space. 

The rest of this paper is organized as follows. In the next section we outline the 
model and calculate the neutrino masses at 1-loop order. Section |3] then presents 
some approximate formulas for neutrino masses and mixing angles, which allow to 
understand how the model can explain the experimental data. We then turn to phe- 
nomenology in section HJ We calculate the decays of the lightest chargino to leptons 
plus scalar neutrino, assuming the (singlet) sneutrinos are the LSP. We compare the 
expected signals with limits on parameters imposed by BR(/i — > e + 7). We then close 
with a short summary. Some formulas for the calculations of loops and LFV decays 
are relegated to the appendix. 



2 Minimal supersymmetric inverse seesaw 
2.1 The model 

The model is defined by the superpotential of the MSSM extended by a pair of singlet 
fields, V c and S with lepton numbers assigned to be —1 and 1, respectively. The total 
superpotential contains then three additional terms [7J 

W = Wmssm + e ab hiL?9 c H b u + M R V C S + ~p s S§ . (2) 
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Note that, in the limit where us — > 0, lepton number is conserved and that the 
parameter Mr does not violate lepton number. We introduce only one generation 
of v c and S. This model is thus the minimal supersymmetric inverse seesaw model 
(MSISM) capable of explaining neutrino data. Previous works used three generations 
of singlets, see e.g. [3|29]. The model conserves R— parity, and as a consequence, the 
lightest SUSY particle is stable. 

With the additional singlet fields the soft SUSY breaking Lagrangian is specified 

by 

-£ soft = -C^ M + mlu^ + mlS^S 

+ [e ah A\l^E h u + B Mr V c S + ±B, S SS + h.c.) , (3) 

where £^| SM contains the usual soft SUSY breaking terms of the MSSM. The param- 
eter B^ s is the analogue of the lepton number violating parameter fig in the superpo- 
tential. The model thus includes two parameters which violate lepton number, both 
will necessarily contribute to the (Majorana) neutrino mass matrix. 



2.2 Tree— level neutrino and sneutrino masses 



From eq. fl2]) we obtain the mass matrix of the neutral fermion fields, which, in the 
basis (u e , u^, u T , u c , S), reads 



M v 















mD 1 





7«d 2 

mo 3 



m Dl m D2 m D[i M R 
M R fl S J 



(4) 



where = h % u v u (i = 1,2,3), with v u being the vacuum expectation value of the 
Higgs field, (H®). For m£,. Mr, one obtains the effective (3 x 3) mass matrix of the 
light neutrinos in the seesaw approximation: 



(M v 

\ m 



mass/*? 



M 



j m Di m Dj 



(5) 
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The lepton number violating parameter us controls the absolute scale of the neutrino 
masses. Eq. (jSJ) shows manifestly the projective nature of the light neutrino mass 
matrix. Thus only one neutrino mass is non-zero at tree-level. However, this result 
is true in general and does not depend on the seesaw approximation. Note also, 
that if rrtOi is of the same order as Mr the correct eigenvalue is found by replacing 
M R ^ M R + J2 m 2 Di in eq. ©. 

The neutrino mass matrix in eq. (jSJ) is diagonalized by an unitary transformation 
in the standard way 

U trT M^U* = dizg(0,0,rn V3 ) . (6) 

In order to obtain a second non-vanishing neutrino mass eigenvalue, loop corrections 
must be included. In this context it is amusing to note that in the non-SUSY case the 
inverse seesaw requires two copies of the singlet fields, v°, Si (i = 1, 2), in order to give 
rise to a viable neutrino mass matrix |16j . even after loop corrections are taken into 
account. 

Assuming CP conservation the 10x10 sneutrino mass matrix can be decomposed 
into two 5x5 matrices for the CP-even, <p R = (V R , V R , V R , V cR , S R ), and CP-odd scalar 
fields, cj) 1 = (?g,£^,z/^, i? c/ ,5 J ), respectively, and reads 




The two mass matrices M\ are given by [7] 

( (M|. + |m|cos2/? + m^.)^ ±(A{ v v u - n m D j cot 0) m D .M R \ 

±{A i hu v u - n m Di cot (3) m 2 uC + M\ + ELi m D k HsM R ± B Mr 

\ m Di M R /j s M r ± B Mr m 2 s + t4 + M R ± J 

where we use a compact form to write these matrices with the index % for the row and 
the index j for the column, i, j = 1,2, 3. The real symmetric mass matrix in eq. (J7]) 
can be diagonalized by a 10x10 orthogonal matrix as follows 

GM 2 G T = di ag (m% i ,...,m%J, (9) 
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with m,~. < . . . < m% . Diagonalizing the mass matrices for the CP-even and CP-odd 

iVi JVio ° ° 

mass matrices M\ separately by 

G ± Ml Gl = diag(m*±, . . . , m?±) , (10) 

leads to a parametrization which is useful for a qualitative discussion of the parameter 
dependence of the neutrino mass matrix which we wish to address below. 



2.3 Neutrino mass matrix at 1-loop order 

We now compute the 1-loop radiative corrections to the neutrino mass matrix. The 
amplitude for the loop contributions to the neutrino self-energy can be generically 
written ajf| 



R 



.fir 



Clearly, the self-energy functions Y7£ v (p 2 ) must be symmetric with interchanging their 
indices due to the Majorana nature of the neutrinos. The 1-loop corrected neutrino 
mass matrix is given by 



M^ 1o °p = m Um {Q)5 mn + Re fe™V) + m Um T^ n {p 



A=0 



(12) 



where m Vm = (0,0, m V3 ) and the self-energy functions Y?gy(p 2 ) are evaluated at p 2 = 
m 2 3 , which is tiny compared to the masses of the particles in the loop, and in excellent 
approximation can be set to zero. The divergences in eq. ffl2l are removed, using the 
minimal subtraction scheme, i.e. by setting the parameter A = 2/(4— <i)— log An = 
0. Here, as usual, d is the number of space-time dimensions, 7^ is the Euler constant, 
and Q is the renormalization scale at which the input parameters are defined. 

The 1-loop improved neutrino mass matrix in eq. ( fl2i) is then diagonalized by an 
unitary matrix denoted as f/ 1-100 ?. The neutrino mixing matrix relating the flavor 
basis (b> a ) and the mass eigenbasis (z/j) of the light neutrinos is then given by 



tr ttI— loop\ 



(U tr U 



(13) 



3 In order to make our results more easily comparable with the case of the standard supersymmetric 
type-I seesaw, we closely follow the notation of [50] . 
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There are two different types of 1-loop diagrams. One class of diagrams exchanges 
Higgses and neutrinos. As we show in detail in the appendix, the flavour structure 
of this loop repeats the flavour structure of the tree-level mass matrix, eq. (151), and 
thus only renormalizes m U3 . More important is the scalar neutrino-antisneutrino loop, 
since it implements a new flavor structure (besides h % u ) and thus generates a second 
non-zero neutrino mass. The new flavor structure is due to the trilinear couplings A\ , 
see eq. ()3]). 

The relevant interaction for the calculation of the self-energy functions is the 
sneutrino-neutralino-neutrino interaction, which is given by the Lagrangian 

C ux o, = xM™jb P R + A L mjb P L )v m N b + h.c. , (14) 



with 



1 

v 

,L 9 



A mjb — --^KUi m N j4 (G b4 - iG b9 ) , (15) 



K jb = --(N* 2 -t a n9 w N; i )(G bl - i G bil+5) )Ur m , (16) 

where g is the SU{2)l gauge coupling and 9w is the weak mixing angle, respectively, 
and N is the unitary 4x4 neutralino mixing matrix, which diagonalizes the neutralino 
mass matrix by N* M x oN~ l = diag(m x o, . . . , rn x o), with m x o > 0. The sneutrino mass 
matrix and diagonalization have been discussed in the previous section. 

The calculation of the self-energy functions then yields 



mn X-j 



~ (4tt) 2 
-1 



^mjb^njb + A m *j b A n J b + A^* b A njb + A^ jb A^* b B (rn x o,m 2 ^J, (17) 



s^mn 
n V2 



(4tt)< 



^rnjb^njb + ^mjb^njb + ^mjb^njb + ^mjb^njb "^1 ( m xJ ' ) ' (^) 



In the limit that the right-chiral couplings A^ b , eq. (fT5l) . are omitted, the result of 
the sneutrino-neutralino loop calculation [30] in the standard type-I (SUSY) seesaw is 
recovered. In the type-I seesaw the couplings A^ b are negligible because of the tiny 
mixing among left-handed and right-handed sneutrino states, which are separated by 
a large mass hierarchy. On the other hand, in the inverse seesaw left-handed and 
right-handed sneutrino states have similar mass scales, such that the contribution of 



S 



couplings A^j b is relevant. As has been noted [30], the functions E™ 2 n are UV finite, as 
the A part in the loop function B^m^, m % ) drops out, because of the orthogonality 
of G. For the same reason is also independent of the renormalization scale Q. 
The off-diagonal functions E™^ n are finite, as expected, and only the 33-element of 
the diagonal elements Ey m gives a non-vanishing contribution to the neutrino mass 
matrix, because of the second term in eq. (Tl2|) . where £y 2 retains a dependence on A 
and Q. Numerically we find that the contributions due to £y 2 , which are multiplied 
by the small neutrino mass, see eq. fll2p . are much smaller than the contributions due 
to E™^ and can be safely neglected in general. The self-energy functions in eq. ( 11 2p is 
then given by the sum of the neutrino-Higgs (see appendix) and sneutrino-neutralino 
contributions 

£-V) = SsiV) + EgSV) , = SviV) + £y 2 V) • (19) 

3 Approximate expressions for neutrino masses and 
fit to experimental data 

The lepton number violating parameters fis and B /Xs govern the scale of neutrino 
physics. B^ s essentially controls the size of the loop contributions, while /is is re- 
stricted due to the tree-level neutrino mass (and thus plays only a sub-leading role in 
the loops). However, only in the limit where both lepton number violating parameters 
vanish, i.e. /is, B fls — > 0, the masses of the CP-even and CP-odd scalars are pair- 
wise equal, i.e. m?+ = m?_, . . . = m?_. In this limit there is then a complete 

1 ^1 ^5 ^5 

cancellation between the contributions of the CP-even and CP-odd scalar loops [T7] . 



The role of other model parameters can be understood with the help of the following 
approximate relations. In the flavour basis one can write the 1-loop contribution to 
the neutrino mass matrix as 

Ml~ loop = a e m e n + b {e m 5 n + 5 m e n ) + c 5 m 5 n , (20) 

with the vectors e m and 5 m defined as 

e m = m Dm M R , (21) 
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A Z V u - V m D m COt (3 , 



(22) 



with a, b and c being coefficients that depend on all other model parameters, see below. 
It is important to note that the e m have the same flavour dependence as the tree-level 
neutrino mass contribution, thus the vectors e m and 5 m must not be aligned in order 
to explain neutrino data correctly. Note also that the structure in (1201) is only strictly 
true, if the soft SUSY breaking parameters are equal for all generations, i.e. 
Mi t = Mi 2 = . Otherwise the new flavour structure introduced by M^. should 
be taken into account. 



The coefficients a, b and c are found with the help of an approximative diagonal- 
ization of the scalar neutrino mass matrices. We give below the formulas for the case 
Mi. 3> Mr (where Mr for simplicity stands for all parameters of the singlet sector) 
since this case is phenomenologically more interesting, as explained in the next sec- 
tion. Formulas for M^. -C Mr can be found easily. Expanding A in the "small" 
parameters e m and 5 m the coefficients read 



E 



m v o 
(4tt) 2 



a j2 

-Jj- ( cos 2 6 + B (m 2 x o,ml+) - cos 2 6L B (m 2 o , m 2 _ ) 
i 

+ sin 2 9 + B Q (m\,m\+) - sin 2 9_B G (m\,m 2 . r 



j 2* 



+st| ,c„^ +B „(^,^)-co^_ Bo K s ,0 



+ sin 2 + Bo(m\,m?~+) - sin 2 6LI? (m 2 , m?- 

Xj "2 X-i "l 



Vj "2 



3* J 

a R a L 



M 2 M R v u 



(sin 29_ ( B (m 2 o , ml- ) - B (m 2 o , ml- ) ) 



+ sin 29+ (B (m x o , m~+ ) - B (m y o , m ~+ ) ) 



(23) 



E 



m v o 



cos 2 +J B o (m 2 o,m 2 
+ sin 2 £ + 5 (m 2 o,m 2 + 

Xj " 2 



cos 2 ^5 (m 2 o,m 2 



X,' 



■■i fsin 2L(B (mJo, m 2 



sin 2 9-B (rrtio, ml- 

Xj v 2 



Bo(m 2 x0 ,ml-)) 
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+ sin 29+ (B (m 2 , m 2 + ) - B (m 2 , m 2 + ) ) 



(24) 



-E 



r w 



jkq ( cos 2 9 + B (m 2 x o, m 2 +) - cos 2 6>_S (m 2 o, m?-) 

L L 



+ sin 2 + B o (mV m 2 +) - sin 2 9_B (m 2 ,m 2 . 



(25) 



Here we have used the abbreviations a J L = —g/2(N* 2 — tan 9 ^N^), g J R = —l/y2Nj4 
and M| = M| + cos 2(3. The mixing angles 9± diagonalize the 2x2 sub-matrices 
of the (v cR , S R ) and iv cl , S^ 7 ) systems, respectively, and are given by 



COS^j 



B 



Mr 



B M~r~ 



[mz 



mi 



sin 9 A 



m 2 c + M\ — m 2 



B Mr 



mi 



,(26) 



fl Bmr . m^ e + — 

C ° S " ^^{m^ - mi, -M\r Sm - ^ H + K_-m^-M^' 1 

where we have neglected the tiny \x$ term in the (2,2) entry and the J2k m D fc term 
(which is of higher order in the seesaw expansion) in the (1,1) entry. The corresponding 
mass eigenvalues are denoted by m 2 + and m?_ . We stress that, in order to derive 

"1,2 u l,2 

the analytic formulas for the sneutrino mixing angles 9± we have implicitly assumed 
that the mixing of the singlet sneutrinos to the left sneutrinos is small, i.e. e m and S m 
are smaller than all other mass squared parameters of the problem. 

Let us consider the case where the mass eigenstates are close to the weak eigen- 
states, i.e. the mixing angles 9± are close to or it/2. This corresponds to the parame- 
ter Bm r being small. It can be shown that the mass squared difference m 2 + —fn 2 - goes 
to zero for cos 9± — > 1 , while for cos 9± — > the mass squared difference of the heavier 
states approaches zero, i.e. m 2 + —m 2 - — > 0. Given this result, from eqs. (1231) - (1251) one 
finds that in the limit Bm r — > 0, only the coefficient a is non-vanishing. For a viable 
neutrino mass matrix, however, we will need also a contribution from the last term 
in eq. (J20l) and this in turn requires therefore a sizeable Bm r - Note that this results 
holds true also for the reversed case, i.e. M^. <C Mr. 
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Figure 1: Two examples of neutrino mass spectra as a function of the parameter B^ s . 
To the left: fi s = 7 eV; to the right = 0.1 eV. 

Neutrino oscillation data require two distinct neutrino mass scales, i.e. the at- 
mospheric and the solar scales. Given the above discussion, in the MSISM neutrino 
masses can be fitted either with one scale generated by tree-level physics, while the 
other is due to the sneutrino-antisneutrino loop or with both scales generated at loop 
level. An example for each case is shown in Fig. [TJ The left panel shows an example 
for the atmospheric scale being due to tree-level physics, with the solar scale generated 
by loops. The right panel shows an example for both masses generated at loop level. 
Which case is realized can not be predicted from the model and depends on the relative 
size of the unknown parameters fis and B^ s . Numerical values used in these figures 
are: M 2 = 500 (GeV), tan/3 = 10, \l = 150 (GeV), M R = M s = 50 (GeV),M„ = 45 
(GeV), m Dl = (GeV), -m D2 = m Da = 2.6 (GeV), M Ll = 700 (GeV), M La = 750 
(GeV), M La = 800 (GeV), B Mr = 50 2 (GeV 2 ), for =e [10~ 2 ,30 2 ] (GeV 2 ). Note 
that on the left panel = 5000 GeV 2 , while in the right panel V5; = 1100 GeV 2 . We 
would like to stress, however, that these are just some random examples. 

Oscillation data fix two Am 2 , namely Am\ TM and AmL but not the absolute scale 
of neutrino masses. Since also the "sign" of Am^ TM is not fixed by oscillation data 
yet, in general three types of spectra can fit solar and atmospheric data. These are 
known in the literature as (a) normal hierarchy; (b) inverse hierarchy and (c) quasi- 
degenerate neutrinos. We note that within the MSISM it is not possible to get all 
three light neutrinos degenerate, thus we will discuss only (a) and (b). 

Any realistic model for neutrino mass must not only explain the absolute values for 
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the atmospheric and solar neutrino mass scales, but also the corresponding leptonic 
mixing angles. As first observed in [3Tj . the so-called tri-bimaximal mixing pattern, 
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(28) 



provides a very good first-order approximation to the measured neutrino angles. This 
pattern can be realized in different ways. However, for normal hierarchical neutrinos 
A^^ iag = (0, m0, Matm) it leads to the following structure of the neutrino mass matrix 
in the flavour basis: 
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(29) 



Here Matm ijn®) represents the atmospheric (solar) mass scale. 

For the inverse hierarchy, _M^ iag = (±Matm,Matm + ^s,0), due to a sign ambi- 
guity in Matm, there are two possible textures 



M HPS 



Matm 



ihi 



(2 





0^ 




(l 


1 l\ 











1 


1 1 


1 


1 


3 


u 


1 


1 J 




^ 1 


1 1 / 



(30) 
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(31) 



Here, ms 



Ami 



-. Comparing eqs. (|29|) -( 131| with the index structure of eq. (120]) 
it is fairly obvious that the MSISM Q can quite easily fit normal hierarchy, whereas 
the case of inverse hierarchy requires a finely tuned cancellation between the different 
contributions (proportional to and 5i) to eq. ( 1201) . We discuss normal hierarchy first. 

Consider the extreme case 6 = 0, see eq. ( l20l) . The structure required by experi- 
mental data could be reproduced with mr, 1 = 0, rri£, 2 = —m£, 3 and W5j = 5 (and vice 



*And, indeed, any model of neutrino mass with this index structure in generation space. 
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versa). However, while the relative importance of the terms a and c can be indepen- 
dently adjusted by adjusting Bm r , b is not independent of a and c at the same time. 
Thus, these equalities are not exact. However, one can use this ansatz as a starting 
point and find valid combinations - within the allowed ranges of neutrino angles - for 
vtld i and Si by a simple iterative procedure. In the numerical scans shown in the next 
section we have always allowed that the range of mj^Jm^. and Si/Sj vary randomly 
within some moderate factor such that all of the allowed range in the neutrino angles 
are covered. With such a random selection of parameters we can fit all angles easily, 
however, there is no prediction and no "typical" size of any neutrino angle. 

The case of inverse hierarchy, however, requires that the tree-level and 1-loop con- 
tribution to the neutrino mass matrix are finely tuned against each other. If, for 
example, we compare the first texture for inverse hierarchy eq. fj30|) with eq. (120]) one 
finds 

ae\ + 2be l 5 1 + cb\ = M ATM (32) 

aeie 2 + b(ei5 2 + ei<5i) + c5i5 2 = — (33) 

o 

2 , r>7 r r2 ^ATM 

ae 2 2 + 2be 2 5 2 + c5 2 = — - — (34) 

i.e. tree-level and 1-loop contributions have to be tuned to cancel each other up to 
J^s j n order to reproduce the desired texture. Similar relations hold for the other 
texture, eq. ( 13"Tj) . We did not attempt to find such fine-tuned solution in the numerical 
scans discussed in the next section. 



4 Lepton flavour violation and collider signals 

In this section we discuss phenomenological aspects of the MSISM. We will concentrate 
on LFV charged lepton decays and the decays of charginos to charged leptons and 
singlet sneutrinos. In general, the new singlets of the MSISM could appear in decay 
chains at the LHC if either (or both) mp or 8 are large, as expected in the MSISM, thus 
potentially altering the phenomenology with respect to MSSM expectations. However, 
the probably most interesting part of the parameter space is that where one of the 
scalar singlets is the LSP, thus being potentially a DM candidate. In this case scalar 
singlets are guaranteed to show up at the end of the supersymmetric decay chains. We 
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will exclusively concentrate on this case in our discussion of chargino decays below. 
Note, however, that LFV £j — > £^ decays are independent of this assumption. 

We consider the decays 

xf^N a + £f, a = l,...,4, ^ = e,/i,r, (35) 

with N\ (N3) being the CP conjugated state to N2 (N4). The relevant piece of the 
Lagrangian for the calculation of the decay widths of (1351) is 

L ix - V = %{C\ ]a P L + C^PR^Na + h.c. , (36) 

with 



c L = -J- 

ija ^ 



gV*(G ai -iG a{i+5) ) -KV*(G a4 -iG aQ ) , (37) 



where the charged lepton Yukawa couplings are = ^ = e,/i,r, and U 

and V are the unitary 2x2 chargino mixing matrices, which diagonalize the chargino 
mass matrix by U* M x ±V~ l = diag(m x ±, m x ±), with m x ± > 0. The decay widths of 
the decays (1351) are finally given as 

mt ^N a + £f) = {m %~ m y (\CL? + \CL?) ■ (38) 

32tt m i v ' 



m 



N 2 



As the members of each CP conjugated pair are always nearly degenerate, 
and ~ rrif fi , they (most likely) cannot be distinguished experimentally. For this 
reason, we sum over the CP-even and associated CP-odd sneutrino states of each CP 
conjugated pair 

r(xt^N 1+2 + ef) = r{xt^Ni + £t) + r(xf^N2 + £f), 

T(xf^N 3+i + £f) = r(xt^N 3 + £t) + T(xt^N 4 + £f). (39) 

To understand the dependence of the decay widths in eq. (J381) on the model parameters, 
one can use an approximate diagonalization of the sneutrino sector as discussed above. 
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If M^, Mr ^> Si,5i, the leading contribution to the decay width to the lightest CP 
conjugated pair, T(xt — > iV 1+2 + if), according to eq. fl38|) is given by 



2 




EIC5 



1L 1 2 
ila I 



(cos 2 9+ + cos 2 6L) , if ^ > 5* 



a=l 
2 




(L 1 2 
ila I 



(cos 2 + cos 2 9-) , if £j <C 5i 



(40) 



a=l 



The results for the decays into the heavier second pair of singlet sneutrino states, 
T(xt — >• N 3+i + if), are obtained by replacing (cos 2 9 + + cos 2 6L) — > (sin 2 9 + + sin 2 9 ) 
in (j40]). 

In our numerical calculations, we have fixed the parameters as follows: M 2 = 700 
GeV, tan/3 = 5, fi = 400 GeV, M u = 700 GeV, M R = M s = M u = 200 GeV, B Mr = 
(200) 2 GeV 2 . This choice is motivated by eq. PU|) which shows that the higgsino 
component of the chargino couples proportional to ef to charged leptons. Other pa- 
rameters have been randomly generated: (I]i m i> i ) 1 ^ 2 £ 10t -4,2 ' 6 ', (Di^ 2 ) 1 ^ 4 £ 10t _4,3 l 
Neutrino data on mixing angles (and mass scales) constrains the other parameters. In 
the numerical examples we adjust the parameters fis and _B MS in such a way that the 
atmospheric neutrino mass scale is determined by the tree-level neutrino mass matrix 
contribution, eq. ([5]), while the solar neutrino mass scale is obtained by the 1-loop 
correction. The component rri£> 1 then has to be considerably smaller than the com- 
ponents rri£) 2 ~ m£>. 3 , so that the reactor neutrino angle is small and the atmospheric 
neutrino mixing angle is maximal; the components Si are all of the same order so that 
the solar mixing angle is large. Note that we have imposed neutrino data to be in 
agreement with the experimental 3a allowed range. Also note that in all the plots we 
have imposed the experimental upper bounds on the low energy LFV radiative decays 



In order to quantify whether the main contribution to the chargino decays is due 
to the parameters mo i or the parameters <5j, in our numerical analysis we define the 
ratio 



We will concentrate on the case where m^. gives the dominant contribution to the 
chargino decay (r > 1). Some comments on the other extreme are given near the end 



br(^- -^e i + 1 ). 




(Ei^ 2 ) 1/4 ' 



(41) 
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of this section. 

For the case r > 1, Fig. [2] shows the correlation of the decay width of the lightest 
chargino to the lightest pair of quasi-degenerate CP conjugated sneutrinos (A^ and 
N 2 ) and a charged lepton £j with respect to the corresponding parameter m 2 D . We 
have checked that this correlation also holds for the chargino decay width T(xt — > 
N 3+ 4 + if), which involves the second lightest pair of quasi-degenerate CP conjugated 
sneutrinos (A3 and N4). This behaviour is as expected from the analytical approx- 
imation in eq. (|40j) . Note, however, that the correlation of the widths involving the 
electron with respect m 2 Di are not as clean than the others. This is due to the con- 
straint on the neutrino reactor angle imposed by neutrino data, which requires rriD 1 
to be much smaller than mo 2 and mo 3 . Comparing the size of these calculated widths 
to typical widths for final states Xi ~ >* X° + W ± and Xi ~ > X° + ^ ±v one finds that 
branching ratios into muon and tau final states can be sizeable, whereas the width to 
final state Aq+2 + e* 2 is expected to be too small to be measurable. 




m Di [GeV 2 ] ra 2 ^ [GeV 2 ] 



Figure 2: Decay width T(xt — > ^1+2 + &t ) of the lightest chargino to the lightest pair 
of quasi-degenerate CP conjugated sneutrinos and a charged lepton = e in the left 
panel and ii = fi in the right panel) as a function of the parameter m 2 Dl (left panel) 
and m 2 D2 (right panel). The plot for m 2 D is very similar to the one for m 2 D2 and thus 
not shown. All plots correspond to the case r > 1, see eq. (|4TI) . 

We note in passing that the product of the decay widths of the lightest chargino to 
one the two lightest pairs of quasi-degenerate CP conjugated sneutrinos and a charged 
lepton ii times the same width but to the charged lepton £j are correlated with the low 
energy LFV process BK(£j — > ici). Again, the correlation involving the electron in the 
final state is less strong than the ones involving only /x and r because of the relative 
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smallness of the parameter imposed by the experimental upper bound on the 
neutrino reactor angle. Since the absolute widths, however, will not be measurable at 
the LHC, more interesting phenomenologically are ratios of partial widths, i.e. ratios 
of branching ratios. 

Fig.[3]shows ratios of branching ratios BR(xf — > N 1+2 + /BR(xf —> ^i +2 + t ± ) 
as a function of BR(/i — > e + 7)/BR(r — > ej) (left panel) and m 2 D2 /m 2 D3 (right panel). 
Again, the same correlations can be found for BR(xf iV 3 +4 + // ± )/BR(xf ^3+4 + 
t ± ). A measurement of both, chargino decays and LFV lepton decays, would therefore 
constitute a consistency check of the scenario we discuss. Note, however, that the 
expected branching ratio for BR(r — > ej) is quite small (at most 10 -12 ) compared to 
current experimental sensitivities. 
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Figure 3: Ratio of branching ratios BR(xf — > N i+2 + /i ± )/BR(xf — > N\ +2 + T± ) as a 
function of BR(/i —7-6 + 7)/BR(r — >• ej) (left panel) and ^n 2 D2 /m 2 Dz (right panel). 



Left panel in Fig. H] shows the correlation of the ratio of branching ratios of the 
lightest chargino decaying to the lightest sneutrino pair and /i divided by its decay 
to the lightest sneutrino pair and r as a function of the atmospheric neutrino mixing 
angle. Recall that these data points have parameters chosen such that the atmospheric 
scale is generated by tree-level physics. The correlation exists for the lightest and for 
the next-to-lightest pair of singlet sneutrinos, if kinematically accessible. Right panel 
in Fig. H] shows the correlation of the ratio of BR(/i — >■ e + 7) divided by BR(r — > e + 7) 
as a function of the atmospheric neutrino mixing angle. As can be seen from both 
panels in Fig. HJ the neutrino sector (the atmospheric mixing angle) is related to 
collider observables (the LFV decays of the lightest chargino to a singlet sneutrino 
and a lepton) as well as low energy LFV observables (the radiative decays of the 
charged leptons). 
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Figure 4: In the left panel, ratio of branching ratios BR(xf — >• iV 1+2 + /i ± )/BR(x5 b — > 
^Vi+2 + r± ) as a function of the atmospheric neutrino mixing angle, tan 2 (#ATM)- In 
the right panel, BR(/i — j- e + 7)/BR(r —> ej) as a function of tan 2 (#ATM)- Both plots 
correspond to the case of r > 1. 

As in our model we have relatively light right handed neutrinos and sneutrinos with 
large Yukawa coupling we should check the contributions from these new particles to 
the muon g — 2. We have calculated the new contributions to a M and verified that our 
numerical points - once they pass the cuts from U — > l 3 1+7 - also pass the experimental 
constraint from a M [35] . 

Finally we would like to comment on the case r 1, i.e. the parameters 5i 
giving the dominant contribution to the neutrino mass matrix and thus to the lightest 
chargino LFV decays. We have scanned the parameter space of the model for such 
solutions and, as expected the decays of the charginos to singlet sneutrinos plus a 
lepton correlate with Si instead of in this extreme. However, in all points we 
have found the absolute widths for the final states N a plus a charged lepton are much 
smaller than in the case r > 1 discussed above (at most of the order of 10~ 5 GeV). 
One expects therefore that the corresponding branching ratios are too small to be 
measured at LHC. 



5 Summary 

The minimal supersymmetric inverse seesaw model (MSISM) with only one pair of 
singlet superfields can explain all existing neutrino oscillation data. We have calculated 
the neutrino mass matrix at 1-loop order and discussed the constraints on model 
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parameters due to the experimentally measured leptonic mixing angles and neutrino 
masses. 

Since in the MSISM one expects the new singlet fields to exist at a mass scale below 
(approximately) TeV, additional phenomenology is expected to show up in experiments 
searching for lepton flavour violation, such as \x — > e^y, and possibly at the LHC. 
Absolute values of branching ratios can not be predicted, but the minimal model 
relates the observed atmospheric angle to some specific ratios of branching ratios. A 
measurement of these ratios can therefore potentially serve as a test of our minimal 
model. 

For the LHC we have concentrated in our discussion on the case that one of the 
singlet scalar fields, a mixture of the scalar neutrino and the scalar singlet S, is the 
lightest supersymmetric particle. This assumption is motivated by the observation 
that this singlet could be the CDM. Charginos can then decay to charged leptons 
plus singlet sneutrinos. A measurement of these decays and low energy lepton flavour 
violating lepton decays, such as fi — > e + 7 and r — > e + 7 would provide an interesting 
test of the minimal supersymmetric inverse seesaw model. 
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Appendices 



A Higgs-heavy neutrino loop 

Here we consider the 1-loop contributions to the neutrino mass matrix which are 
mediated through the Higgs-heavy neutrino loops. For their calculation, the required 
Lagrangian is given by 

C HOuuL = -hiE* {k+m Ut^P L v m H° u + h.c. , (A.l) 

where E denotes the unitary mixing matrix, which diagonalizes the mass matrix of the 
neutral fermion fields of eq. (TjJ, E*M U E~ 1 = diag(m7vj = diag(0, 0, m„ 3 , m Ll , m L , 2 ). 
We note that only the contribution of the two heavy neutrinos is significant, and 
thus in ( 1A.1I) uj^, k = 1,2, denotes the two heavy neutrinos which are a mixture 
of the fermionic states v c and S. Contributions of light neutrinos in the loop are 
negligible, since the corresponding mixing elements E^, i = 1,2,3, are tiny. The 
neutral component of the up-type Higgs doublet can be expressed in terms of mass 
eigenstates as follows [32] 

Hi = v u + —= [cos ah + sin aH° + z(cos f3A° + sin f3G )} , (A.2) 
V2 

where h°, H° are the two CP-even scalar fields, with m h o < m H o and the corresponding 
mixing angle a, A is the CP-odd scalar field, while G° is the Goldstone field, with 
tan/3 = v u /vd- From the calculation of the self-energy functions of these contributions 
we obtain 



£51" = ^^^^>^™^(f + 3)45oK fc ,m^) , (A.3) 
S ™ = ^^^^^^r^l^ + 3)4| 2 i?iK fe ,m^) , (A.4) 

where we have introduced the shorthand notations % r = (cos a, sin a, i cos (3, i sin /3) 
and rrifi r = (m h o,m H o,m A o,mz), to sum up the various Higgs boson contributions. 
In eqs. (lA.3[) - f1A.4[) . the standard 2-point loop integrals B (x,y) and Bi(x,y), when 
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evaluated at zero momentum (p 2 = 0), can be written as 



B (x,y) = A + l + log 

. . 1 
B 1 (x,y) = -- 



Q' 2 



o 2 

A + l + log — 



X X 

log- 

x-y y 

x 



X 

log- 



V x-y y 



(A.5) 
(A.6) 



We see from eqs. flA.3p -( ]A~~4]) that the flavor structure is determined by the product 
of the neutrino Yukawa couplings h™h2 (in the flavor basis), and has therefore the 
same structure as the tree-level contribution, see eq. (J5]). As a result, only the 33- 
element in E™" and receive non-vanishing contributions. Thus, by including only 
the Higgs-heavy neutrino loop contributions still two of the light neutrinos remain 
massless. 



B LFV lepton decays 

Here we summarize the formulas for the calculation of the two-body LFV lepton decay 
rates in the MSISM with only one generation of singlet superfields. The formulas are 
derived from the superpotential in eq. (]2]) and the soft SUSY breaking Lagrangian in 
eq. (E]). In the context of the SUSY inverse seesaw mechanism with three generation 
of singlet superfields and mSugra boundary conditions see Ref. 



The gauge invariant amplitudes of the decays £j (p) — > l i (p — q) + 7(g), lj 
/i, r; £j = e, fi, can be defined as 



T = z'ee*(g)w^.(p - q)[a^ 'q v {a L ^P L + a R ^P R )}u tj {p) . (B.l) 

In the calculation of the left and right amplitudes, ul.r-, we neglected terms propor- 
tional to the small lepton mass m^.. The heavy lepton contributions give rise to the 
right amplitude as 



hl i9 2 ^ P * w / -4g_+45g_-33g + 10 3s 3 



with s = m 2 Nk /m 2 v . Note that the loop function in eq. ( 1B.2I) coincides with eq. (68) 
of [53], and differs by a constant, —5/6, from eq. (B.2) of Ref. [31]. For £~(p) — > 
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^iip — o) + 7(9) decays this constant does not contibute, due to the unitarity of 
the coupling matrices, but for g — 2 the main contribution for light neutrinos comes 
precisely from this constant and therefore the correct loop function for both cases is 
eq. ( 1B.2|) . The sneutrino-chargino loop contributions to the right amplitude read 



sc 

T R,ij 



10 



167T 



EE 

fe=l a=l 



^ika^jka 2 

Na 



m- 



t 2 - 5t - 2 
,12(t-l) 3 

/ t-3 



tint 



2(t-l) 4 
hit \ 



Hka^jka 2 



2(f-l) 2 



(B.3) 



with £ = m^±/m| and the couplings defined in eq. ( |37j) . The sneutrino-chargino loop 



contributions to the left amplitude are obtained from the right ones by interchanging 
the left and right chiral couplings, i.e. crf^- = a R c i:j (L <H- R). With the definition of 
the amplitude in eq. ( IB. II) . the corresponding decay widths are given by 



a 



SC \2 
'L,ij\ 



+ 14% 4 



HL |2\ 



(B.4) 



with a = e 2 /47r and where again in the kinematics we have neglected the terms 
proportional to the small lepton mass m^. 



C Muon anomalous magnetic moment 

The formulas for the muon anomalous magnetic moment can be derived from those in 
the previous section. Defining, as usual, 

we have from eq. fIB.lj) 

a M = -2m^ (cr L>2 2 + 0X22) (C.2) 

where &l,22 an d ctr,22 are defined, for our model, in eq. (IB.2p and eq. (IB.3[) . The 
sum in eq. (1B.2[) should run only over the heavy neutral leptons, as the light neutrino 
contribution is alreday included in the Standard Model. We have checked that our 
formulas reproduced the PDG [35J value for the sum of the W and Z diagrams, at 
1-loop order in the Standard Model. 
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